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Abstract — The performance of free-space optical (FSO) com- 
munication systems is compromised by atmospheric fading and 
pointing errors. The pointing errors are widely considered as a 
combination of two components: boresight and jitter. A statistical 
model is investigated for pointing errors with nonzero boresight 
by taking into account the laser beamwidth, detector aperture 
size, and jitter variance. A novel closed-form probability density 
function (PDF) is derived for this new nonzero boresight pointing 
error model. Furthermore, we obtain closed-form PDF for 
the composite lognormal turbulence channels and finite series 
approximate PDF for the composite Gamma-Gamma turbulence 
channels, which is suitable for terrestrial FSO applications 
impaired by building sway. We conduct error rate analysis of 
on-off keying signaling with intensity modulation and direct de- 
tection over the lognormal and Gamma-Gamma fading channels. 
Asymptotic error rate analysis and outage probability of such a 
system are also presented based on the derived composite PDFs. 
It is shown that the boresight can only affect the coding gain, 
while the diversity order is determined by the atmospheric fading 
effect as well as the pointing error effect. 

I. Introduction 

In free-space optical (FSO) communication systems, laser 
beams are transmitted through atmosphere, ranging from hun- 
dreds of meters to several kilometers in terrestrial solutions 
and to thousands of kilometers in earth-to-satellite links [1]. 
Such systems are vulnerable to both atmospheric attenuation 
and pointing errors, which can severely deteriorate its error 
rate performance. 

In terrestrial FSO communication systems, the transceivers 
are often positioned at the top of tall buildings to obtain a line 
of sight. In such systems, the atmospheric turbulence, building 
sway, building vibration and thermal expansion of building 
can degrade the link performance. The atmospheric turbulence 
can cause random temporal and spatial irradiance fluctuation 
in the optical beam. The building sway, building vibration 
and thermal expansion of building can result in pointing 
errors that consist of two components: boresight and jitter. 
The boresight is the fixed displacement between beam center 
and center of the detector. Although typical terrestrial FSO 
systems are initially installed with near zero boresight error, 
the boresight is still considerable due to the thermal expansion 
of the building. In [2], the pointing errors are recorded over 
three days in a TerraLink laser communication system. The 
radial displacement shows a cyclical pattern every 24 hours, 
which suggests the thermal expansion of the building plays an 
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important role in pointing errors. The reported boresight can 
be as high as 0.3mrad [2]. In [3], the values of boresight in 
the range of 0.0039 — 0.0117mrad are used. In [4], the author 
considers boresight up to 0.03mrad in a horizontal FSO link. 
The jitter is the random offset of the beam center at detector 
plane, which is mainly caused by building sway and building 
vibration. The typical value of jitter standard deviation in a 
terrestrial FSO system is below 0.3mrad [2]. In [3], the jitter 
standard deviations of 0.0033 — O.OlOOmrad are used. In [4], 
the author considers jitter standard deviation up to O.Olmrad. 

In satellite-to-ground and intersatellite communications, the 
transmitter and receiver have high relative velocity, and there 
is mechanical noise due to satellite-based motion and gimbal 
friction [5]. Thus, it is difficult to realize perfect tracking, 
while jitter and boresight can also arise as residual pointing 
error. 

Numerous studies on pointing error effects have been re- 
ported. In [6], the author proposes a mathematical model 
to minimize the transmitter power and the beam divergence 
angle in an urban optical wireless communication system with 
pointing errors caused by building sway. In [7], a maximum- 
likelihood estimator is developed to estimate the boresight and 
jitter component of the pointing error. In this system model, a 
point detector and nonzero boresight component are assumed. 
In a follow-up work [4], the same authors further consider the 
effects of atmospheric turbulence for lognormal and Gamma- 
Gamma fading channels, and they adopt a wave-optics based 
approach to evaluate the channel capacity. In [8], a statistical 
pointing error model is proposed by considering the laser 
beamwidth, jitter variance, and detector size. In this model, a 
closed-form probability density function (PDF) of the pointing 
error loss factor is derived by assuming that the pointing error 
has zero boresight. In [8], the size of detector area is taken 
into consideration. When the size of the detector area becomes 
negligibly small, this pointing error loss factor coincides with 
the one derived for the point detector in earlier works [6], 
[9]. Furthermore, the authors derive the optimal beamwidth 
to maximize the channel capacity at a given probability of 
outage. In [10], the error rate performance of FSO links 
over if-distributed turbulence channels in the presence of 
pointing error is studied by assuming intensity modulation 
direct detection (IM/DD) with on-off keying (OOK). The 
system under study uses the same pointing error model de- 
veloped in [8]. The average bit-error rate (BER) is presented 
in closed-form using the Meijer G-function. In a related work 
[11], the authors study the BER performance of a heterodyne 
differential phase shift keying optical wireless communication 
system in the presence of pointing error over the Gamma- 
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Gamma turbulence channels. In [12], two optimization models 
are proposed to mitigate the pointing error effects with zero 
boresight by taking into account the transmitter power, wave- 
length, transmitter and receiver telescope gains. In [13], Al- 
ary pulse-position modulation is investigated with impairments 
from atmospheric turbulence and pointing error with zero 
boresight. In [14], asymptotic BER performance is analyzed 
for FSO communication systems using transmit laser selection 
over atmospheric turbulence channels with the same pointing 
error model developed in [8]. A statistical channel model is 
recently proposed for multiple-input multiple-output (MIMO) 
FSO communication systems over atmospheric fading chan- 
nels with pointing errors [15]. Both outage probability and 
diversity order are studied, and it is found that the diversity 
order is determined by pointing error effects other than the 
number of transmitters or receivers. More recently, the point- 
ing error model in [8] is generalized by modeling the radial 
displacement distance with a Hoyt or Nakagami-g distribution 
[16], which allows the horizontal and elevation jitters to have 
nonidentical jitter standard deviations. 

The pointing error model developed in [8] is widely used in 
the literature [10], [11], [13], [14]. In this model, the boresight 
component of pointing error is assumed to be zero, and both 
horizontal and elevation displacements are assumed to follow 
an independent, identically distributed zero-mean Gaussian 
distribution [9]. As a result, the random radial displacement at 
the receiver is Rayleigh distributed. In our work, we assume 
nonzero boresight component and assess the performance of 
terrestrial FSO systems. We derive a statistical model to 
describe the pointing error effects with nonzero boresight 
components by considering a finite size circular detector. We 
obtain a novel closed-form PDF of our generalized pointing 
error model. Our pointing error PDF specializes to the one in 
[8] when the boresight error is zero. Furthermore, we derive 
a closed-form PDF for the composite lognormal channels 
and an approximate series PDF for the composite Gamma- 
Gamma channels using our new pointing error model. Our 
series PDF for the composite Gamma-Gamma channels is 
applicable when 7 2 > a, where 7 is the ratio between the 
equivalent beamwidth and jitter standard deviation, and a 
is a parameter related to the large scale turbulence cells in 
Gamma-Gamma turbulence fading model. We find that the 
condition j 2 > a is satisfied in many terrestrial FSO systems 
impaired by building sway. Finally, we study both the error rate 
performance, asymptotic error rate performance, and outage 
probability for an OOK IM/DD based FSO system in the 
presence of nonzero boresight pointing errors. 

The remainder of this paper is organized as follows. Section 
II describes our system model and derives the closed-form 
PDF of our pointing error model. We review two important 
atmospheric turbulence channels and derive closed-form PDF 
for the composite lognormal channels and series PDF for 
composite Gamma-Gamma turbulence channels when -f 2 > a. 
We present several examples where j 2 > a is satisfied 
and show that -f 2 > a is valid in many terrestrial FSO 
systems impaired by building sway. In Section III, we derive 
expressions of BER, asymptotic BER, and outage probability 
of an FSO system over the lognormal and Gamma-Gamma 



turbulence channels with nonzero boresight pointing errors. 
Section IV presents some numerical results, and Section V 
makes several important conclusions. 

II. Turbulence Channels with Nonzero Boresight 
Pointing Errors 

OOK modulation is widely used in practical FSO systems 
[24], [25] since optical communication systems with higher 
order modulation are complex to implement. In this work, we 
consider an OOK FSO system with IM/DD. The transmitter 
modulates the data bits directly onto the intensity of an optical 
beam, which propagates through an atmospheric turbulence 
channel in the presence of pointing errors. Since the channel 
coherence time is on the order of msec and the data rate is 
assumed to be on the order of Gbps, we can therefore adopt 
a slow fading channel model. The received optical power is 
converted into an electrical signal through direct detection at 
the photodetector. Assuming additive white Gaussian noise 
(AWGN) for the thermal/shot noise and unit detector respon- 
sitivity, we can express the received signal y at the detector 
as 

y = hx + n (1) 

where x is the transmit intensity being either or 2P t , 
where P t is the average transmitted optical power, h is the 
channel gain, n is zero-mean AWGN with variance a 2 n . Fol- 
lowing [8], [11], [15], we assume that the off-axis scintillation 
changes slowly near the spot of boresight displacement, and 
use a constant value of scintillation index to characterize 
the atmospheric turbulence. Therefore, the pointing error and 
atmospheric fading are independent, and the channel gain h 
can be expressed as h = hih p h a , where hi represents the 
path loss which is a constant at given weather condition and 
link distance, h p is the pointing error loss factor, and h a is 
the atmospheric fading loss factor. Note that pointing error 
loss factor h p and atmospheric fading loss factor h a are both 
random variables (RVs). 

A. Nonzero Boresight Pointing Error Model 

When a Gaussian beam propagates through distance z from 
the transmitter to a circular detector with aperture radius a, 
and the instantaneous radial displacement between the beam 
centroid and the detector center is r, the fraction of the 
collected power at receiver can be approximated as [8] 



h p (r; z) w A exp 



2r 2 



it' , 



(2) 



where A is the fraction of the collected power at r = 0, and 
w Zcq is the equivalent beamwidth. We have A = [crf(w)] 2 
and w 2 = w 2 .yZ_ r /^.2\ . where v = ^/tt/2— is the 



z 2v cxp( — v 2 ) ' 

ratio between aperture radius and beamwidth, and erf (2) = 
^= f£ e - * dt is the error function. The beamwidth w z can be 
approximated by w z = 6z, where 9 is the transmit divergence 
angle describing the increase in beam radius with distance 
from the transmitter. For example, at a range of 1 km, a 1 
mrad divergence produces a beam radius of lm at receiver. 
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The approximation in (2) is accurate when w z /a > 6, which 
is satisfied in typical terrestrial FSO communication systems. 

At the receiver aperture plane, we can express the radial 
displacement vector as r = [r x ,r y ] T , where r x and r y , 
respectively, denote the displacements located along the hor- 
izontal and elevation axes at the detector plane. We consider 
a nonzero boresight error in addition to the random jitters, 
and model r x and r y as nonzero mean Gaussian distributed 
RVs, i.e., r x ~ Af(/j, x ,a x ), r y ~ M{p y ,cr 2 ). Then the radial 

displacement r = |r| = yV 2 + r 2 follows the Beckmann 
distribution [26] 

fr(r) = 



2^a x u y 



r-2-K 

/ exp 
Jo 



(r cos <p — ji x ) (r s\n (f> — jiyY 



2al 



2a 2 , 



d<j>. 



(3) 



The Beckmann distribution is also known as lognormal-Rician 
distribution [27], [28], which is used to describe the PDF of 
fading channels in general. It is a versatile model which applies 
to a variety of distributions. For examples, it can specialize to 
Rayleigh when fi x = [i y = 0, a x = a y ; Rician with ^i 2 +fj, 2 ^ 
0, a x = G y \ Hoyt distribution when \i x = \x y = 0, a x ^ 
o y [16]; and single-sided Gaussian when [i x = p, y = a x = 
0,(7^0 [29]. 

In satellite FSO communication systems, it is widely ac- 
cepted that the jitter variance is the same for both horizontal 
and elevation axes [21], [22]. In terrestrial FSO systems, 
however, the jitter is mainly caused by turbulence and building 
motion. Since the turbulence cells randomly appear on the 
beam path, and the building might be considered to sway in 
orthogonal and parallel directions to the beam path with equal 
probabilities, we can therefore assume a 2 = a 2 = a 2 [8], 
[11], [15]. As a result, the PDF of radial displacement r in 
(3) becomes Rician 



fr(r) 



exp 



jr 2 + s 2 ) 
2a 2 



Io 



(4) 



where s = [i 2 + ^ 2 is the boresight displacement, and / o (0 
is the modified Bessel function of the first kind with order zero. 
From (2) and (4), we derive the PDF of nonzero boresight 
pointing error as 



fhjhp) = 



7 2 exp 



(3) 



A 1 




In 



< h p < A 
(5) 

where 7 = w Zcq /2a s is the ratio between the equivalent 
beamwidth and jitter standard deviation, which is a measure 
of the severity of the pointing error effect. If we consider 
zero boresight error with s = 0, our pointing error model in 
(5) specializes to the one in [8]. Our analytical result in (5) 
is accurate when w z /a > 6. In Appendix A, we compare our 
derived PDF of h p in (5) with the exact PDF, which is obtained 
numerically, without using (2) to show the accuracy of our 



analytical model. In Appendix B, we derive the moments of 



hp as 



E [ K ] = 



A n "i 2 



exp 



n + 7 2 (™ + 7 2 )2ct 2 

where E[] denotes the expectation. 



(6) 



B. Atmospheric Fading Channels 

In order to study the effects of turbulence-induced fading 
on the performance of FSO systems, several statistical models 
characterizing atmospheric fading have been proposed. For 
weak turbulence conditions, we use the lognormal fading 
model to characterize the atmospheric fading h a whose PDF 
is given by 



fh a (h a ) 



1 



2h a y // 2na 



exp 



(liih a + 2a 2 x Y 
8a% 



(7) 



where a\ is the log-amplitude variance given by a\ k> 
ax/4 = 0.3lk 7/6 C 2 z 11/6 [8], where a 2 R is the Rytov variance 
for a plane wave, C 2 is the index of refraction structure 
parameter of atmosphere, and k = 2tt/X is the optical 
wavenumber with A being the wavelength. The parameters of 
the lognormal fading model can be measured directly for FSO 
systems [30]. For medium to strong turbulence conditions, we 
use the Gamma-Gamma turbulence model to characterize the 
atmospheric fading h a whose PDF is given by 



fh a (h a ) = 



2(a/3) 



(«+/3)/2 



~h a 



(2^pfCj 



(8) 



where T(-) is the Gamma function, and X ct _ ( g(-) is the 
modified Bessel function of the second kind of order a — /3. 
The parameters a and (3 are related to the small scale and 
large scale eddies respectively. 

Depending on the value of the Rytov variance, we can 
approximately categorize the turbulence regime as follows 
[32]: the weak turbulence regime (o-fj < 0.3) and the moderate 
to strong turbulence regime (o\ > 0.3). 

C. Composite PDF with Generalized Pointing Error 

The PDF of channel gain h = hih p h a can be calculated as 
[8] 



fh(h) = / fh\h a {h\h a )f ha (h a )dh a 



~I h a hi fh '{ 



h a hi 



(9) 



fh a (h a )dh a . 



In weak turbulence conditions, we use lognormal fading 
model to characterize the atmospheric turbulence fading. We 
substitute (7) into (9), and after some mathematical manipula- 
tion, the composite PDF of lognormal fading with the pointing 
error model in (5) can be obtained as 



2(A a h i r \ u c ) 
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where u a = ^ + 2o 2 x -f + 2a 2 x j 4 , u b 



z eq 

4ct| 7 2 , u c = ^8 +<Tx)> and erfc ( u ) = 1 - crf ( u ) 

is the complementary error function. It is worthy to note that 
our PDF in (10) can specialize to eq. (14) in [8] when we set 
the boresight to zero (s = 0). 

Since hi is a constant, while h p and h a are independent 
RVs, The nth moment of h can be obtained as E[h n ] = 
h?E[h%]E[h2}. With the derived moments for h p in (6) and 
the moments of the lognormal RV, we have 



E[h n ] 



n + 7 2 

x cxp ( -2a x n + 2a x n 2 



ns 

(n + 7 2 )2<7 s V ' 

(11) 



For medium to strong turbulence conditions, we model h a 
as a Gamma-Gamma distributed RV. After some mathematical 
derivation shown in Appendix C, we obtain a finite series 
approximation of the composite PDF as 



f GG (h) « f GG (h) = 
1 { ap^ j 




(vj{a,P)hP- 1+i - Vj(P,a)h a - 1+i ) 



(12a) 



where 

Vj(a,P) 



r 



exp 



o 2 2/2 

' 2<j1 2/3~2j 2 +2j 



) sin" 1 ((a- /3)tt) 



T(a)T(P)T(j - (a - 0) + 1)| - (J3 - 7 2 + j) 



and J = L7 2 — aj- Therefore (12a) is an approximation 1 of the 
composite PDF of Gamma-Gamma fading with pointing error. 
Our composite PDF in (12a) is applicable under the condition 
7 2 > a. Using a rough approximation, we can show that the 
condition ~/ 2 > a corresponds to ^§ > max{6, 2.31a r 5 }- 
Now we present some examples where this condition holds: 
When crfj = 5.0 and the jitter angle is 0.3mrad, the transmit 
divergence is larger than 1.28mrad. When u\ = 25.0 and the 
transmit beam divergence is 2mrad, the jitter angle is less than 
0.35mrad. In terrestrial FSO systems with link distance less 
than 5km, the system parameters have the following typical 
values: 5 ~ 20cm for the receiver diameter, 2 ~ lOmrad for 
the transmit divergence, and <~ 0.3mrad for the boresight 
and jitter angle [2], [17]. It can be shown that such systems 
generally operate under the condition 7 2 > a. However, in 
applications with long link range or narrow beam divergence, 
our PDF in (12a) is not applicable since -f 2 > a may not 
hold in such applications. In the rest of paper, we assume 
7 2 > a unless otherwise stated. In (12a), the parameters a, (3 
are required to satisfy (a — /3) ^ Z. When (a — ft) e Z, 
one can add a small value e to a to satisfy the condition 

'This finite series approximation is accurate for strong turbulence con- 
ditions, but it can be inaccurate for weak turbulence conditions. A large 
boresight and jitter can also make the series inaccurate. 



(a — /3) ^ Z [33]. The value of e is empirically chosen, which 
should be small (i.e. 10~ 3 ) that will not change the Gamma- 
Gamma distribution dramatically. For example, if we have a = 
2.12,13 = 1.12, we can add 0.001 to a, and use a = 2.121 in 
the Gamma-Gamma model. 

Similar to the moments obtained for composite lognormal 
channel, using (6) and the moments for the Gamma-Gamma 
RV we can obtain the nth moment of h as 



E[h ri 



(A /i;)"7 2 r(q + n)r(/3 + n) 
(n + 7 2 )r(a)r(/?)(a / 3) r 



exp 



ns' 



(n + 1 2 )2a\ 



(13) 



III. Error Rate Performance 



Following the system model described in Section II, the 
BER of IM/DD with OOK modulation conditioned on the 
channel gain is 



Pe(e\h) = Q 



Pth 



— erf c I — — 
2 W2cr„ 



(14) 



where Q(x) = exp ^— du is the Gaussian Q- 

function. The average BER can be obtained as 



f 

Jo 



P e (e\h)f h (h)dh. 



(15) 



According to (1), we have E[x 2 ] = 2P 2 . Assuming a slow 
fading channel, the instantaneous signal-to-noise ratio (SNR) 



is defined as SNR = 



2Pfh 2 



. We define the average SNR as 



(12b) t ^ le rece i ve d SNR when there is no fading and pointing error 

V ' ,rn[,.2l ,> ™ , ,»rn ^ cTivrrT 2P? 



(E[h 2 } = 1). Thus, we have the average SNR as SNR 



A. Bit-Error Rate 

1) Composite lognormal model: We derive the average 
BER of lognormal fading with nonzero boresight pointing 
error by substituting (10) into (15). The BER expression can 
be written as 



eXN — 



7 2 exp(n a ) 



x^^-Wc( ln ^ + Mb Vrfcf 



Pt 



\V2a n 



h dh. 



(16) 



Using a change of variable rule, eq. (16) can be expressed as 



D TV I 2 \ 

Pe,LN = — CXp (U a - 7 li b ) 



I 



x / exp( 7 u c x)erfc(x)erfc 



P t A hi 



\[2o n exp(uft — u c x) 



dx. 
(17) 



By introducing an auxiliary parameter B and partitioning the 
integration interval in (17) into [— 00, B] and [-B,oo], we can 
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e.LN 



l 2 u c 



exp (u a 



l 2 u b 



) x / cxp (7 

J — 00 



2 u c a;)erfc (x) crfc ( ^*^ ^' cxp (u c x — b) 



\f2a n 



dx + Rb 



(18) 



rewrite (17) as (18), where R B is the approximation error 
given by 



Rb = -^7^ cxp (u a - j 2 u b ) 



x / exp(7 2 u c a;)erfc(a;)erfc 



V2a n exp(u b - u c x) 



(19) 



In Appendix D, it is shown that the approximation error Rb 
can be upper bounded by 



Rb < 

y/TT'J 2 1 



exp [u a - 7 u b + 



A 1 

IK 



erfc B 



7 2 u c 



We quantify Rb under various system parameters and calcu- 
late the values of B for different Rb- The results are presented 
in Appendix D, showing that the approximation error Rb 
decreases rapidly with increasing B. We can always adjust 
the value of B to make R B arbitrarily small. Therefore, eq. 
(18) can be accurately approximated as 



P e ,LN ~ P e ,LN 
,-B 



l 2 u c 



cxp (u a - 7 2 u b ) 



X 



/exp(7 2 u c x)erfc(x)erfc I —= — 
-00 *" \\/2(T r , 



PtAoh 



exp(u b — u c x) 



dx. 
(21) 



Using a series expansion of the complementary error function 
[34, eq. (06.27.06.0002.01)] 



2 ~ (-l) fc *2*+i 
erfc(z) = 1 - _ ^ 



7T fc!(2fc + l) 



and an integral identity [34, eq. (06.27.21.0011.01)] 

M„~t„{„~\ a~ _ ^ ( J>z 



I 



e bz erk(az)dz — ^ (e bz eric(az) — el^erf az 



(23) 

we derive an infinite series expression of the BER in (21) as 
(24) at the top of next page, where 

S3 = exp(( 7 2 + 2j + l)u c B)crk(B) 

+ exp(( 7 2 + 2j + l)\ 2 /4)crfc(( 7 2 + 2j + l)u c /2 - B). 

(25) 

The infinite series in (24) can be rigorously shown to be 
convergent, and a detailed proof is presented in Appendix E. 

2) Composite Gamma-Gamma Model: For the Gamma- 
Gamma fading with nonzero boresight pointing error, we 
substitute (12a) into (15) to find its BER. Using an integral 
formula [34, eq. (06.27.21.0132.01)] 



J 



1 /a + 1 

/ttq; V 2 



(26) 



we can obtain the approximate BER in terms of a finite 
series as (27). We comment (27) is obtained without additional 
approximations over those in (12a). The approximation error 
defined as e(P e ) = \P e —P e \ is discussed in Appendix F. When 
P t is beyond the minimum required value in Table IV, we can 



dx guarantee the relative error v D e ' ; and 



e(Po„t, gg) 



less than 



P C ,GG Pout.GC 

10 -6 , where P ou t,GG denotes the outage probability of system 
in the composite Gamma-Gamma fading case and e(P outi GG) 
denotes the approximation error \P ou t.GG — Pout.Gcl- With the 
setting s/a = 1.0 and a s /a = 1.0, the minimum required P t 
is below —10 dBm, which can be satisfied for practical FSO 
systems. 



(20) B- Asymptotic Error Rate Analysis 



At asymptotically high SNR, average symbol error probabil- 
ity of an uncoded system in fading channels can be accurately 
approximated as P^° = (G c ■ SNR) _Gd , where Gd is the di- 



versity order indicating how fast the BER decreases with SNR 
in high SNR region, and G c is the coding gain determining 
the shift of the BER curve in SNR relative to the benchmark 
curve (SNR)~ Gd . The asymptotic BER P£° can reveal the 
behavior of BER in high SNR region, which is helpful in 
conceptual understanding of performance limiting factors in 
communications over fading channels [36]. The diversity order 
and coding gain are determined from instantaneous SNR's 
PDF through its behavior near the origin. The PDF of channel 
gain can be expanded into power series as [36] 



lim if h (h) =ah t +g t (h) 

/i-s-0 



(28) 



(22) where g t (h) satisfies lim S*M = Q. Considering SNR 



2P (hi 

— ^2 , the diversity order and coding gain are obtained as 
[36f 

G d = (29) 



and 



G c = 



-2VoT(| + l)' 



(30) 



In the following, we first derive the power series expansion 
of the PDF near its origin, then we obtain the diversity order 
and coding gain from (29) and (30). The asymptotic BER 
is obtained from P^° = (G c ■ SNR)~ Gd . For the composite 
lognormal fading with pointing error, we obtain the PDF of the 
channel gain h near the origin by applying the limit operation 
to (10) directly, the result is shown as 



lim f L N(h) 

h— >0 



-exp(u a )/i 7 1 +g J 2_ 1 (h). 



(31) 



IEEE TRANSACTIONS ON COMMUNICATIONS 



6 



Pe.LN — 



l 2 u c 



exp (u a - l 2 u b ) 



l 2 u c 



exp('j 2 u c B)erk(B) + exp 



4 2 

7 X 



erfc 



7 V 



J=0 



\ 2j+1 exp(- Ufc (2 J + l)) 



j!(2j + l)( v y2a/°N (2j + l+ 7 > c ^ 



(24) 



e.GG 



Pe.GG = 



j=0 



2 a r(2^|±i) 

a+j 



13 + J 

(IP 2 
Vi(J},a) -f 

\ 71 



(27) 



From (31), we obtain the diversity order as Gd — \ — ~^§r- 
This indicates that the diversity order is determined by the ratio 
between equivalent beamwidth and jitter standard deviation. 
More specifically, the boresight component of pointing error 
does not affect the diversity order in lognormal fading channel. 
The coding gain can also be obtained from (31) as 



G c = 



2^ 2 " 1 r(^ + i)cxpK) 



(32) 



With SNR 
rate as 



2Pf 



, we can present the asymptotic bit-error 



poo 

r e,LN — 



1?- X Y ( £ + i) exp (u„) / 2 p2 



V (33) 

For the Gamma-Gamma composite fading channels, we 
derive the power series expansion of PDF near the origin in 
Appendix C, and it is given by 

lim f GG {h) = v (a, ^h - 1 + gp-i (h) (34) 

where v (a, (3) follows the definition in (12b) and we have 
assumed 7 2 > a. From (34), we obtain the diversity order as 
Gd = (3/2 and the coding gain as 



2 /»-l r 



exp 



2<J? 



2£-2 7 2 



a[3 
A a hi 



Y{ a )T{P) sin((a - P)n)T(-(a - /?) + 1)| 7 2 - P\0 



x r 



(2±i) ^ 



(35) 



Thus, the asymptotic bit-error rate can be expressed as 

poo 

^e,GG — 



r(a)r(/3) sin((a - /3)7r)r(-(a - 0) + 1)| 7 2 - /3|/3 ( 36 > 



Note that the diversity order is Gd = (3/2 when 7 2 > a, 
which implies that the Gamma-Gamma fading effect is more 
dominant than the pointing error effect with respect to the 
BER performance at high average SNR region. However, when 
7 2 < a, the diversity order will depend on both 7 2 and 
the boresight s. An explicit expression of diversity order for 
nonzero boresight pointing case is difficult to obtain. For the 
special case with zero boresight, the diversity order is given 
as min{ 7 2 /2,/3/2} [37]. 

To evaluate the performance loss caused by boresight error, 
we define an SNR penalty factor in dB as 



SNR boresight = 10 log 10 



SNRpoc 

,nonzcroborcsight 

SNRf 



.zcroborcsight 



(37) 



which represents the constant SNR loss caused by boresight 
error at certain error probability when SNR is asymptot- 
ically large. From (33) and (36), we obtain the SNR penalty 
factor for the composite lognormal and the composite Gamma- 
Gamma fading, respectively, as 



SNR, 



boresight, LN — 



80 
In 10 



and 



SNR boreisi9 ; lt , GG 



40 
In 10 



Aa 2 (3 



(38) 



(39) 



where we have assumed 7 2 > a, that is w 2 > 4er 2 /3 
in (39). The results in (38) and (39) show that, for both 
composite lognormal and composite Gamma-Gamma cases, 
the boresight error imposes a constant SNR penalty on the 
error rate performance when SNR is large. As expected, 
larger beamwidth oo Zeq can mitigate the adverse impact of the 
boresight error. 



C. Outage Probability 

In a slow fading FSO system with IM/DD OOK, we 
assume the receiver has perfect knowledge of channel gain, 
and the data rate of transmitter is R bits/channel use. For 
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Fig. 1. BER performance of IM/DD OOK over the lognormal fading with pig. 2. BER for the composite lognormal channel (a\ = 0.01, s/a = 2) 



zero and nonzero boresight pointing errors. 



with different jitter values. 



OOK signaling, the instantaneous capacity corresponding to a 
specific channel realization h! is given by [8] 

c(h') = jYl fy\*(y\x)p x {x) iog 2 fj j^y- d y < 4 °) 

where x e {0, 2P t }, p x (x = 0) = p x (x = 
2P t ) = 0.5, f y \ x (y\x) = 7 =L=exp((^#^), f y (y) = 

^2 x p x ( x )fy\x(y\ x )- We define the outage probability as the 
probability that the instantaneous capacity C(h) is not suffi- 
cient to support the data rate Rq, and it is given by [8] 

P mt =Prob(/i<C- 1 ( J R )) (41) 
which can be calculated as 

Pout = I ° fh(h)dh (42) 
Jo 

where h a = C~ 1 (R ). Therefore, we can use the composite 
PDF, which has been derived for different turbulence models, 
to calculate the outage probability at a given transmission 
rate. For weak turbulence condition, substituting (10) into (42) 
and using the integral identity in (23) we obtain the outage 
probability as (43) at the top of next page. For medium to 
strong turbulence conditions, by substituting (12a) into (42), 
we obtain the outage probability in (44). 

It is worthy to mention that we can obtain the diversity 
order using the outage probability derived in (43) and (44), the 
derivation follows that in [15]. The results are G^ln = 7 2 /2 
and Gd,GG = (3/2, which coincide with the ones obtained 
using the power series expansion approach. 

IV. Numerical Results 

In this section, we adopt the system settings shown in 
Table I, which are used in many practical terrestrial FSO 
communication systems [2], [17], [39]. Under two typical 
weather conditions shown in Table II [8], we carry out the 
error rate performance of an FSO communication system. In 
order to investigate the degradation effects induced by the 
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Fig. 3. BER for the composite lognormal channel (cr^ = 0.01, s/a = 
2, u s I a = 1.5) with different beamwidth values. 
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Fig. 4. Outage probability of an FSO system over the lognormal fading with 
zero and nonzero boresight pointing errors. 
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out,LN 



7 f ln irt 

exp(u a )erfc 



hp 
A h t 



+ exp u a - u & 7 



u. 



In h ° 



+ u b 



(43) 



Pout.GG ~ Pout.GG = 

exp 



r( a )r(/3)sin(( a -/3)7r)^ o 



E 



a/? 
Aoh t 



( 



2f3-2~f2+2j 



L0+3 



T(j - (a + - (/? - 7 2 + 3W + J) ° 



2 2/2 

2q-27 2 +2j 



,a+j 



r(j + (a - 0) + l)j!| - (a - 7 2 + j)\(a + j) 



(44) 



TABLE II 

Sample system parameters under light fog and clear sky weather conditions 



Weather 


Visibility 


Path loss (hi) 


CI level 


Fading model 


°l 


Parameters 


Light fog 


0.5 km 


0.008 


10 -15 m -2/3 


Lognormal 


0.05 
0.2 


a\ = 0.0125 
a\ = 0.05 


Clear sky 


10 km 


0.9 


10 -14 m -2/3 


Gamma-Gamma 


0.6 
2.0 


a = 5.41,/? = 3.78 
a = 3.99,13 = 1.70 




Fig. 6. Outage probability of an FSO system over the Gamma-Gamma fading 
Fig. 5. BER performance of IM/DD OOK over the Gamma-Gamma fading with zero and nonzero boresight pointing errors, 
with zero and nonzero boresight pointing errors. 



TABLE I 
System Settings 



Parameter 


Value 


Receiver Diameter (2a) 


20 cm 


Noise standard deviation (<r n ) 


cr„ = 10^ 7 A/Hz 


Link distance 


1 km 


Transmit Divergence at 1/e 


1 mrad 


Corresponding beam radius (u> z ) 


100 cm 


Jitter angle 


0.1 mrad 


Corresponding jitter standard deviation (<r s ) 


10 cm 



boresight error, we use two different values of normalized 
boresight displacement s/a = 0,3. 

In weak turbulence regime, we study the BER of an FSO 
link over the lognormal fading channel with nonzero boresight 
pointing errors. We calculate the exact BER from (16), the 
approximate BER from (24), and the asymptotic BER from 
(33). The BER curves are plotted in Fig. 1 against the 
transmitted optical power P t . From Fig. 1, we can see how the 
boresight displacement s affects the BER performance of the 
system. We can infer from the result that for a Gaussian beam 
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Normalized boresight displacement h p x1(T 3 



Fig. 7. The SNR penalty factor induced by boresight error in different Fig. 8. Comparison of the analytical PDF in (5) and the exact PDF of h p 
turbulence conditions. under various system settings (w z /a = 10). 



the energy collected at receiver aperture decreases with in- 
creasing boresight errors, and therefore the BER performance 
deteriorates with a large boresight displacement. However, if 
we assume a tophat beam profile, the boresight would not 
affect the collected power at the receiver as long as the 
boresight displacement s is smaller than the difference of 
distance between the beamwidth at receiver and the receiver 
aperture. From the asymptotic BER curves in Fig. 2 and Fig. 3, 
we find that the diversity order is determined by the equivalent 
beamwidth as well as the jitter variance. A similar finding is 
reported in [15] for a MIMO FSO system with zero boresight. 
With larger beamwidth or smaller jitter variance, the diversity 
order of the system becomes larger. The outage probability of 
an FSO system with code rate Ro = 0.5 (bits per channel use) 
over the composite lognormal fading channel is presented in 
Fig. 4, where we use (43) to calculate outage probability. It 
can be seen that the outage probability of the system worsens 
with increasing boresight errors. 

In medium to strong turbulence regimes, we study the 
BER of an FSO link over the Gamma-Gamma fading channel 
with nonzero boresight pointing errors. We calculate the 
exact BER from (15), the approximate BER from (27) (with 
J = L7 2 - aj), and the asymptotic BER from (36). The BER 
curves are plotted in Fig. 5 against the transmitted power P t . 
It can be seen that our series solution developed in (27) can 
accurately approximate the exact BER when P t is beyond a 
certain threshold (see Appendix F). However, for small values 
of P t , the series approximation in (27) can be inaccurate. This 
is the limitation of the series approach for the Gamma-Gamma 
channels with pointing errors, and the same limitation can also 
be seen from Fig. 2 of [16]. As expected, the BER performance 
worsens when the boresight displacement s becomes larger. 
Moreover, the boresight error causes a horizontal shift of 
the BER curve, resulting in an SNR penalty factor for the 
error rate performance. From the asymptotic curves in Fig. 5 
we can find that the boresight displacement does not affect 
the diversity order of the system. Assuming a code rate of 
Ro = 0.5 (bits per channel use), we present the exact outage 



probability obtained from (C.l) and (42), and the approximate 
outage probability calculated using (44) in Fig. 6, and it can be 
seen that the approximate outage probability is accurate. (The 
threshold of P t for accurate outage probability approximation 
is shown in Appendix F.) In Fig. 7, we plot the SNR penalty 
factor versus the normalized boresight displacement s in both 
composite lognormal and composite Gamma-Gamma cases. 
The result indicates that the boresight displacement has larger 
penalty factor on SNR when the turbulence is weaker. 

V. Conclusions 

A nonzero boresight pointing error model was investigated 
for a Gaussian laser beam propagating through atmospheric 
channels by taking into account the beamwidth, detector 
aperture, and jitter variance. Our derivation is based on the 
assumption that the boresight component of pointing error 
effects is not negligible. A closed-form PDF was derived 
for the nonzero boresight pointing error model. We derived 
closed-form composite PDF for the lognormal fading. We 
also derived an approximate PDF for the Gamma-Gamma 
fading under the condition -y 2 > a, which is suitable for 
terrestrial FSO systems impaired by building sway. Highly 
accurate convergent series BER expression was obtained for 
the lognormal fading. The series approximation of BER over 
the Gamma-Gamma channels can lead to inaccurate estimation 
if the transmitter power is not sufficiently large. Based on 
the asymptotic analysis, we observed that the boresight error 
causes an SNR penalty factor on error rate performance at high 
SNR. By examining the asymptotic BER curves, we found 
that the diversity order of the FSO system over the composite 
lognormal fading channel is solely determined by the pointing 
error parameter -f 2 and the boresight component does not 
affect the diversity order. While in the composite Gamma- 
Gamma fading channel, the diversity order is determined by 
either the Gamma-Gamma fading effect or the pointing error 
effect. 
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Appendix A 

Analytical and Numerical Results of f hp (h p ) 

When a Gaussian beam propagates through distance z from 
the transmitter to a circular detector with aperture radius a, 
the fraction of the collected power at receiver is [8] 



h p {r;z) = G(r) 

a p\f a 2 —x' 2 



a J — y/a 2 — x 



2-/2 TTWt 



■oxp( -2 ( *' r) l + yl2 )dy'dx' 



wt 



(A.l) 



where r is the instantaneous radial displacement between the 
beam centroid and the detector center. We use the approxima- 
tion given in (2) to derive the analytical PDF of h p in (5), and 
the approximation in (2) has normalized mean-squared error 
less than 10~ 3 when w z /a > 6. Using (A.l), we obtain the 
exact PDF of h v as 



Since there is no explicit expression of G~ 1 (h p ), we calculate 
it numerically. In Fig. 8, we compare the analytical PDF of 
hp in (5) as well as the exact PDF of h p without using (2) 
under various system settings. It is shown that our analytical 
model in (5) is accurate even with large boresight (s/a = 5) 
and large jitter (a s /a = 3). 



Appendix B 
Derivation of the moments of h v 



The moments of the generalized pointing error h p are given 



by 



E[h;] = [ Ao h;.f hp (h p )dh p . (B.i) 

Jo 



Substituting (5) into (B.I), we obtain 



E 



v - I dhp 



(B.2) 

where v = \j — ^ — — . Using a series representation of 
Jo(-) [35, Eq.(8.445)], we have 

E [ K ] 

7 2 exp(-^) ^ unW _^ 1 ( SV \^ 



Jo p hwfWv 



dhp. 



(B.3) 



Since each term of the series is non-negative and the infinite 
series uniformly converges to Iq (§|)> we can swap the 



integral and infinite summation and write (B.3) as 

E [ K ] 

7 2 cxp(-f o 



A 7 



00 i r-A a 

(-4) 



m=0 

7 2 exp 



h nW-l (^.Z^hA\ dh 



AT 



E 

m=0 



8<rf 



(m\y 



Jo 



dx 



(B.4) 



Applying an integral identity [35, Eq.(4.294.10)] to (B.4), we 
have 



h p (hp) = f r (G-\h p )). dG ^. (A.2) E [K\ 



7 exp 



A ol 



2a? 



exp - 



m=0 

2 

ns 



(n + 7 2 )2a 2 



( i 2 * 2 V 



1 



(n + 7 2 ) 



(B.5) 



Appendix C 
Gamma-Gamma Composite PDF 

1 ) Composite PDF Approximation: We substitute (8) into 
(9) and write the composite PDF of the Gamma-Gamma 
fading with nonzero boresight pointing error as (C.l) at 
the top of next page. App lying a change of variable rule 

x = ^1 y^- In Aai^ )' ec l- (C.l) can be expressed as 

(C.2), where W(h) = br^^K^p ^y^exp (^)). 
Using a series expansion of the modified Bessel function of 
the second kind [34, eq. (03.04.06.0002.01)] 



K v {x) = 



7T 



oo r 

y 

2sin(7rw) ^ 



(x/2) 



2p—v 



(x/2) 



2p+v 



_T(p-v + l)p\ T(p + v + l)p\ 

(C3) 

where it requires v £ Z and \x\ < oo, we can express W(h) 
as 



W(h) = 



T(j - (a - 0) + 



2sin(7r(a-^)) 



E 

3=0 



r(j + (a-p) + m 



-h? 



) 



+Q-1 



(C4) 
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= 2 T *„ P (-»v f ;)M)<^>/V 
(Mh) 1 r(o)rG9) 



/ /la 1 
Jh/(A h,) 




■In 



ft 



haAoh^ 



K a -f3(2^af3h a )dh a 



(CI) 



8 7 2 exp(- S 2 /2a, 2 )(a/3) (a+/3)/2 
(A ^)^r( a )r(/3) W 2 



Jo 



.TCXp 



2x 2 fa + f3 



-7 



4 7 r 7 2 exp(- S 2 /2 f x 2 )M)(» + ^ 2 fa(a+ff)/ 
(A o / li ) (a+0)/2 r(a)r(/3) sin((a - P)M . 



2-1 



E 

J=0 



a/3h 

r(j -(n- I) + \ \AofH 

1 / afih 

'T(j + (a-/3) + l)j\ \Aohi 



f 

JO 



XeX P ( - ^ 2 + i) ) ^0 i^X ) d.r 



j+(a-/3)/2 -oo 



/ 



X cxp — 



(a - 7 2 +i) ] /:: ( S ..r ) (' 



(C5) 



Substituting (C.4) into (C.2), and after some manipulation, we 
have (C.5). 

In the following derivation, we use an integral identity [35, 
eq. (6.643.2)] 



Jo 



x u -^e- ax I 2v (2/3^)dx 

E£±£i> !«-»m.„,. 

T(2v + 1) \2a ) \ a 



(C6) 



where M u/V (-) is the Whittaker function. Hence, by using 
another the identity [34, eq. (07.44.03.0041.01)] 



Mr, 



i t rn{z) = exp z 1 ™ m\ I cxp(z) - ^2 



m— 1 i. ' 

Z K 



k=0 



k\ j ' 

TO > 

(C7) 

we can represent (C.5) as (12). In our derivation, the integral 
identity (C.3) requires that (a — f$) d_ Z, and the integral iden- 
tity (C.6) requires that the summation index j in (C.5) must 
be lower than or equal to J = L7 2 — ot\ , where [a;J denotes 
the largest integer not greater than x. Such truncation method 
is also used in [16] to estimate the BER of systems over the 
Gamma-Gamma channels with Hoyt distributed pointing error. 
Due to the truncation of the infinite series and the constraint on 
the number of terms j < [j 2 — a\ , our series PDF in (12) may 
not converge to the exact PDF for weak turbulence. However, 
beyond certain threshold of transmit power P t (shown in 
Appendix F), the approximate BER is accurate since the PDF 
focih) near the origin can be accurately described by (12). 

2) PDF Near the Origin: To obtain the power series 
expansion of the Gamma-Gamma composite PDF near the 



origin, we can express W(h) as 
lim W(h) = 

h— >0 

2sin(7r(a-/3))r(-(a-/3)-f 



(C8) 



1) 



From (C.8) and (C.2), we can derive the PDF near the origin 

as 



lim f G a{h) = 

h—>0 



2f3-2-{ 2 



)sin- 1 ((a-/3) 7 r)^_ i 



r( a )r(/3)r(-( a -/3) + i)|-(/3- 7 2 )| 

+ 5/3-1 W- 



(C9) 



Appendix D 
Bound on approximation error 

By applying an upper bound erfc(a;) < exp(— a; 2 ) to (19), 
we can upper bound Rb as 



Rb < 

7V 



cxp (u a - -f 2 u b ) 



00 



p2 

,2„. „ _2 r t a2t,2 . 



x / cxp I j"u c x — x" — — ± 1 rA Q h l exp(u c x — b) \dx 

JB V 2<T n / 

< - Mc exp (u a — 7 2 m&) / cxp {^ 2 u c x — x 2 )dx. 
4 J B 



(D.l) 
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Using an integral identity [35, Eq.(2.33.1)] 

exp(— (ax 2 + 2bx + c))dx = 



1 HF (b 2 - 
-W-exp 

2 V a \ a 



erf V aa; H 1= 



(D.2) 



we derive the upper bound for R B as 

Rb < 

V^1 2 U C ( 2 / 

exp \u a - 7 M h 



2 2 

7 



erfc B 



(D.3) 

For different Rb, the required values of B, which are cal- 
culated by (D.3), are shown in Table III, where we set 
w z /a — 10, s/a = 1.0, a s /a = 1.0. It is found that the 
approximation error Rb decreases rapidly with increasing B, 
and therefore we can adjust the value of B to make R B 
arbitrarily small. 

Appendix E 
Proof of convergence of series 

The series in (24) contains the expression 

Ps = 



2 

7^ 



! ^ {-iy ( p t Aphi 

%t^' ! (2j + 1) V V2a n 



2j+l 



exp(-t/. b (2j + l)) 
(2j + l + 7 2 K 

(E.l) 



We now use the ratio test to assess the convergence of this 
series. The absolute ratio between two consecutive terms is 



a j+1 



& + !) (%7) 2 exp(-2 Wfe )(2j + 1 + 7 2 ) Sj+1 



(E.2) 



(. 7 + l)(2. 7 +3)(2. 7 +3 + 7 2 ) Sj 
By applying the bounds of the erfc(-) function [41] 
1 



1 - — ^ ) cxp(-a; 2 ) < erfc(x) < — L- — ■' J 
2x 2 / Jttx 



exp(— x ) 
(E.3) 



to Sj+i (upper bound) and Sj (lower bound), we have 



< 



b + (7 2 +2j+ 1 3)i^-s exp(2i?u c ) 



B 2B 2 ) + 



( 7 2 +2j + l)^-B 



2(( 7 2 +2j + l)^-B) 2 / 

(E.4) 

Applying the limit operation to both sides of (E.4) and noting 
that Sj and Sj+i are positive, we have 

1 



< lim < 



l 

2B 2 



fij 1 

For the rest of the terms in (E.2), we have 



(E.5) 



{2j + 1} exp(-2 U6 )(2j + 1 + 7 2 ) 
lim 7T tt^4 ^vT^rr ^ = 0. 



0- + 1)(2j+3)(2j+3 + 7 2 ) 



(E.6) 



Since limits in (E.5) and (E.6) both exist, we have 

dj+i 



lim 

j-s-oo 



lim 



(2j + l)(^) 2 (2i + l + 7 2 ) 



• lim 



^+1 



j^oocxp(2u fc )(j + l)(2j + 3)(2fc + 3 + 7 2 ) j^oo Sj 
= 0. 

(E.7) 



Thus, we have lim c ^ ±1 
Ps in (24) is convergent. 



0. We conclude that the series 



Appendix F 
Approximation error 

After the numerical evaluation of the approximation error 
in (27) under different system parameters, we find our series 
solutions in (27) and (44) may not approach to the exact results 
when P t is small. This is a limitation of the series approach. 
In Table IV, we show the minimum required P t that can 
guarantee the relative error < 10~ 6 and £( p°" t) < 10~ 6 

e "out 

for three representative turbulence conditions. From Table IV, 
we observe that the minimum required P t becomes larger with 
increasing boresight or jitter. 
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